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ABSTRACT 



Context. The Sun and other cool stars harbouring outer convection zones manifest magnetic activity in their atmospheres. The con- 
nection between this activity and the properties of a deep-seated dynamo generating the magnetic flux is not well understood. 
>A/Vt?s. By employing physical models, we study the spatial and temporal characteristics of the observable surface field for various 
stellar parameters. 

Methods. We combine models for magnetic flux generation, buoyancy instability, and transport, which encompass the entire convec- 
tion zone. The model components are: (/) a thin-layer aQ. dynamo at the base of the convection zone; (//) buoyancy instabilities and 
the rise of flux tubes through the convection zone in 3D, which provides a physically consistent determination of emergence latitudes 
and tilt angles; and (///) horizontal flux transport at the surface. 

Results. For solar- type stars and rotation periods longer than about 10 days, the latitudinal dynamo waves generated by the deep- 
seated aQ dynamo are faithfully reflected by the surface distribution of magnetic flux. For rotation periods of the order of two days, 
however, Coriolis acceleration of rising flux loops leads to surface flux emergence at much higher latitudes than the dynamo waves at 
the bottom of the convection zone reach. A similar result is found for a KOV star with a rotation period of two days. In the case of a 
rapidly rotating Kl subgiant, overlapping dynamo waves lead to noisy activity cycles and mixed-polarity fields at high latitudes. 
Conclusions. The combined model reproduces the basic observed features of the solar cycle. The diff'erences between the latitude 
distributions of the magnetic field at the bottom of the convection zone and the emerging surface flux grow with increasing rotation 
rate and convection zone depth, becoming quite substantial for rapidly rotating dwarfs and subgiants. The dynamical evolution of 
buoyantly rising magnetic flux should be considered as an essential ingredient in stellar dynamo models. 

Key words. Sun: activity — Sun: dynamo — stars: interiors — stars: late- type — stars: activity — stars: magnetic field 



1. Introduction 

It is commonly accepted that the large-scale magnetic field of 
the Sun is maintained by a hydromagnetic dynamo process, 
which operates in the convective envelope and leads to an oscil- 
latory magnetic field with an average period of 1 1 years. Among 
the processes considered in the current the oretical context for 
. the magnetic activity cycle (se e reviews by IOssendrijverll2003l: 
' ICharbonneaul 120051: iFanI l2009l) are: (1) generation of toroidal 
magnetic field from poloidal field by dififerential rotation (Q- 
eflTect), (2) formation and instabilities of toroidal magnetic flux 
' tubes, (3) buoyant rise of toroidal magnetic flux, (4) generation 
of poloidal from toroidal magnetic field by helical convection or 
through twisting of rising flux loops by the Coriolis force (a- 
eflTect), and (5) transport of magnetic flux by meridional flow. 

In rapidly rotating stars with outer convection zones, 
manifestations of magnetic activity are observed throughout 
the electromagnetic spectrum. Distribution and coverage of 
starspots can be quite dififerent from the solar activity patterns 
dStrassmeier 2009): spots near the rotational poles, and high UV 
and X-ray fluxes indicative of a large surface coverage of mag- 
netic regions. Observations of stellar magnet ic activity provide 
constraints of stellar dynamo mode ls (see e.g. JStrassmeieil20Q5l: 
ISchussledl2005l: lBerdvuginall2005h . Hence, it is necessary to in- 
vestigate the links between magnetic flux generation and trans- 
port, as well as the eflTects of stellar structure and rotation on the 
spatio-temporal distribution of the surface field. 



It is often tacitly assumed that the latitudinal distribution 
of toroidal magnetic fields generated in stellar interiors can be 
taken to represent the surface emergence patterns. However, 
the correspondence between the latitudinal distribution of the 
dynamo-generated field and the emergence pattern is not self- 
evident, because (a) the stability properties of flux tubes at 
the bottom of the convection zone depend sensitively on field 
strength and their position in latitude (Ferriz-Mas & Schii ssleJ 
fT993rfT995h . and (b) rising flux tubes in a rotating star are de- 
flected towards the poles, owing to angular momentum conser- 
vation (IChoudhuri & GilmanI [19871 ISchiissler & Solankilll992l: 
iGranzer et al.ll2000b . 

Numerical simulati ons of the rise of flux tubes in the 
convection zone (e-g.. iD'Silva & Choudhuril 119931: iFan et al.l 
119941: ICaligarietal]ll995Lll998h have successfullv reproduced 
many observed properties of sunspot groups, among which the 
tilt angle (the angle between the line connecting the oppo- 
site polarity regions and the local latitudinal circle) is of par- 
ticular importance for Babco ck-Leighton-type dynamos (e.g., 
iDikpati & Charbonneaulll999l) . Here, we consider the connec- 
tion between the dynamo mechanism operating in the stellar in- 
terior, the transport of toroidal magnetic flux through the con- 
vection zone, and the emerged flux evolving under the efifects of 
near-surface flows. We combine models for three processes: (1) 
a dynamo operating in the overshoot layer at the bottom of the 
convection zone, (2) the magnetic buoyancy instability and rise 
of magnetic flux tubes through the convection zone, (3) the trans- 
port of magnetic flux at the surface. In this first exploratory study 
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we have used a very simple one-dimensional dynamo model that 
yields the time-latitude characteristics of the solar cycle. A more 
complete two-dimensional dynamo model will be employed in a 
forthcoming study. 

Prelimin ary results obtained with our combined model were 
presented by llsik etaP (l20Q 7a). In this paper, we give a detailed 
account of the model setup, discuss the effects of different rota- 
tion rates in Sun-like stars, and apply the model to a rapidly ro- 
tating KO-type main sequence star, and a Kl-type subgiant star. 



2. Dynamo model 

We consider the generation of magnetic flux in the overshoot 
layer at the bottom of the convection zone (at radius ro - 
0.73/? q). We assume a kinematic afl dynamo operating in a thin 
layer (ISchmitt & Schiisslerll 19891) . The model is based on the as- 
sumption that the radial diffusion of magnetic flux out of the 
layer is partly compensated by downward flux pumping by con- 
vective flows. This thin-layer aQ. dynamo exhibits magnetic cy- 
cles and latitudinal propagation of dynamo waves based upon 
a radial gradient of the rotational angular velocity. Differential 
rotation in latitude and meridional circulation, although proba- 
bly relevant for the operation of solar/stellar dynamos cannot be 
consistently incorporated into this simple model. We chose to 
use this model as a simple means of obtaining solar-like time- 
latitude behaviour of the toroidal field at the bottom of the con- 
vection zone and to avoid overloading this exploratory study 
with the complications (and uncertain parameters) of a more 
complex model. These models can of course be incorporated into 
our model framework, as we intend to do in future studies. 



the a-effect. For the Sun, we take ro = 5.07 x 10^^ cm, which 
corresponds to the middle of the convective overshoot layer (ac- 
cording to the stratification model used in Sect. [3]), and kro = 3, 
so that a quarter of the wavelength of the radial dependence of 
the magnetic field corresponds roughly to the thickness of the 
convection zone. The second term on the r.h.s. of Eq. ^ rep- 
resents the buoyant loss of toroidal magnetic flux from the dy- 
namo layer. The quantity is the critical mean field strength, 
above which the flux is lost from the layer with a characteristic 
timescale r, corresponding to the growth time of the magnetic 
buoyancy instability of toroidal flux tubes (Sect. 13.21) . The value 
of ^0 scales with the critical field strength for the instability (av- 
erage value between 0° and 40° latitude), which depends on the 
given stellar model and rotation rate (see Figs. [5^, [7^, [9^, [TOh , 
[T2h). The nonlinear function / describes the ^-dependence of 
the flux loss in the form 



/l|-l = 



-sgn(B)'(B/Bo-lf 




ifB>Bo 
otherwise. 



(5) 



We rewrite the dynamo equations in nondimensional form by 
taking ro as our unit of length, the diffusion time, r^/rj, as the 
unit of time, ^o as the unit field strength, ao as the maximum 
absolute magnitude of the Of-effect, and ro^o as the unit of the 
vector potential. Using the same symbols as before, (A, B), for 
the nondimensional quantities, Eqs. (O and © transform into 



dB 
dt 



Rn — (Asme)^Q-f(B) 



1 



. . dB^ 

sm^ — - 
sin6>^6'\ de 



sin^6/ 



(6) 



2.1. Dynamo equations 

We use spherical polar coordinates (r, 0, 0) and decompose the 
azimuthally averaged magnetic field, B, into a toroidal compo- 
nent, (0, 0, B(r, 6, 0), and a poloidal component described by the 
vector potential, A = (0, 0, A(r, 6, t)), yielding 



B 



B{r, e, t)e^ -h V X (A(r, 6, t)e^) , 



(1) 



— = RaB+ - — — 
dt sinOde 



sm^— - — — 
do [sin^e 



(7) 



The induction effects are characterised by two dimensionless 
numbers, namely a Reynolds number for the a-effect 



Rrr — 



aoro 



(8) 



where is the unit vector in the azimuthal direction. Following 
iHoyng et"aD (Il994l) . the radial dependence of A and B is as- 
sumed to have the form of a spherical wave 



A(r,0,t) = A(0J)— Qxpiikr), 
r 

B(r,e,t) = B(e,t)—exp(ikr), 
r 



(2) 



reflecting flux loss from the thin dynamo layer at r = ro by mag- 
netic diffusion in the radial direction. The mean magnetic field at 
r = ro as a function of the colatitude, 6, and time, t, is governed 
by the dynamo equations for A{6, t) and B(6, t), viz. 



f =Q'w|(Asin^?) + ^/(|- 
dt do T \Bo 
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r 



1 d 



sinOde 



. dB\ B 
sm 6— — 

do sin^e 



- (krofB 



dA 7] 
— = a{e)B+^ 
dt rl 



1 d 



sin do 



. dA\ 



sin^6> 



- (krofA 



(3) 



(4) 



where r/ is the turbulent magnetic diffusivity, ^\6) is the radial 
gradient of the angular velocity at r = ro, and a(6) represents 



which represents the efficacy of induction by the Of-effect, and a 
Reynolds number for the radial shear (Q-effect) 



Ra = 



(9) 



which represents the efficacy of induction by the differential 
rotation. The dimensionless number 2, which determines the 
strength of the flux-loss nonlinearity is given by the ratio of the 
timescales for diffusion and flux loss 



T 7/T 



(10) 



In the dynamo model, we consider only one hemisphere and 
assume dipolar parity for the magnetic field, i.e., antisymmetry 
with respect to equatorial plane. The boundary conditions at the 
pole and at the equator are thus given by 



5 = A = 0for6/ = 

^ dAsmO ^ ^ 
5 = 0, — = for 6/ = 71 2. 

de 



(11) 

(12) 
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We set the turbulent difFusivity to 7/ ^ 2.96 x 10^^ cm^ s"^. For 
the solar- type model, ao ^ -11.7 cm s"^ leads to a magnetic 
cycle period of 22 years. 

Equations ©-((Tl) with the boundary conditions Eqs. (fTTT) - 
(fT2l) are nu merically solved with an i mplicit finite-difference 
method (cf. ISchmitt & Schusslerlll989l) using 90 grid points in 
latitude. 



2.2. Differential rotation and a-effect 

The radial rotational shear in the dynamo layer, fl\6) = 
da(r, 6)/dr\r =ro, is taken in the form 



Q\e) = -%pI(cos 0) = --^(3 cos^ -h 1), 



(13) 



where P2(cos0) is the associated Legendre polynomial. 
Equation ([T3]) with = 10"^^ cm"^ s"^ roughly represents the 
shear pr ofile in the solar ta chocline as determined by helioseis- 
mology (ISchou et al.lll998h . 

The latitudinal dependence of the Of-eff'ect is assumed as 



a(0) 



ao sin [n (6 - Oq) / (n/2 - Oq)] : for ^ ^ 







: for6><^o, 



(14) 



where ao is the amplitude, and = 55°, which corresponds to 
the zero crossing of the function Q\6). The functions a(6) and 
fl\6), normalised to their amplitudes, are shown in Fig.[T] In or- 
der for the dynamo waves to propagate equatorward, one must 
have afl' < 0. To fulfil this criterion, ao has been chosen to be 
negative, because the sign of radial shear in the solar tachocline 
is positive for > 55°. The choice of Eq. ([T4l) and the sign of 
ao are mot ivated by the a-effect d riven by magnetic flux tube in- 
stabilities dFerriz-Mas et al and unstable magnetostrophic 
waves (ISchinittii2QQ3l) . 



3. Emergence of magnetic flux 

We assume that (1) the toroidal magnetic field in the dynamo 
layer undergoes a magnetic Rayleigh-Taylor instability, which 
leads to the formation of magnetic flux tubes (iFanit2QQU ,2QQ9). 
(2) the flux tubes reac h a mechanical equilibrium s tate by devel- 
oping an internal flow oreno-Insertis et al Ii992|),and(3)their 
equilibrium location is at the middle of the overshoot region. The 
flux tubes become subject to the undulatory (Parker) instability 
once their field strength exceeds a critical value, Bcr, which is a 
function of latitude. We determine the lin ear stability properties 
following iFerriz-Mas & Schiisslerl (Il995l) and carry out numer- 
ical simulations of the nonlinear evolution of the unstable flux 
tubes and their rise towards the surface. 



3.1. MHD equations for magnetic flux tubes 

The equation of motion for the material inside a flux tube, in a 
reference frame r otating with the angular velocity of the tube, O, 
can be written as (iFerriz-Mas & Schiisslerll 19931) 



B^\^(B-V)B 



An 



+pi^g - n X (O X r)j -h 2piVi X n -h /^D, 



(15) 



where DjDt = d/dt-\-V'V is the Lagrangian derivative. The sub- 
script / denotes quantities inside the flux tube and, in the follow- 
ing, the subscript e denotes external quantities. The terms on the 



right hand side of the equation are, respectively, the total pres- 
sure force (gas and magnetic), magnetic tension force, effective 
gravity (including the centrifugal force), Coriolis force, and hy- 
drodynamic drag force. For the drag force, we use the expression 
for a flow past a straight spherical cylinder. 



PeV±V±_ 



(16) 



where i;^ = is the perpendicular component of the relative 
velocity of the tube with respect to the external medium, is 
the unit vector along v±, Cd is the hydrod ynamic drag coefl&cient 
which we take as unity (iBatchelodI 19671) . pe is the density of the 
external medium, and Rt is the cross-sectional radius of the tube. 

The combination of the equations for continuity and mag- 
netic induction in ideal MHD leads to Walen's equation 



(17) 



For the energy equation, we assume that the mass elements of 
the flux tube evolve isentropically. This is justified because the 
timescale for heat exchange by radiation in the deep convec- 
tion zone is much long er than the timescale of flux tube motions 
Im oreno-lnsertisi 19861) . The system of equations is closed by the 
equation of state for an ideal gas and the condition of magnetic 
flux conservation. 

We consider the equations gi ven above in the framework of 
th e thin flux tube approximation (|SDruitlll98lb . in the form given 
by lFerriz-Mas & Schussled(ll993r 



3.2. The escape of magnetic flux tubes from the dynamo 
layer 

Magnetic flux tubes leave the dynamo layer owing to the undula- 
tory buoyancy instability, which sets in for B > B^r and leads to 
flux loops rising through the convection zone. To determine the 
number, field strength, and latitude of unstable flux tubes, we 
assign a probability, p, for a flux loop to start rising at a given 
latitude and time. 

In the case of the Sun, we h ave B^y - 10^ G in the middle of 
the convective overshoot layer (ISchiissler et al.l [19941) . The cor- 
responding magnetic energy density is about 100 times higher 
than the kinetic energy density of convective motions. We as- 
sume that the shear-generated mean toroidal field, B, is frag- 
mented by the magnetic Rayleigh-Taylor instability into flux 
tubes of strength Byy, such that 



(18) 



where / is the filling factor of flux tubes contributing to the mean 
field. For Byi = 10^ G and / = 0.1, we have B ^ 10^ G, which is 
of the order of the equipartition field strength, Bq^ = {Anpy^Vc, 
where Vc is the convective velocity determined by using a mix- 
ing length model. The individual flux tubes are considered to be 
intensified to field strengths of the order of 10^ G by a combina- 
tion of rotational shear and the conversion of potential energy in 
the external stratification as suggested by iRempel & SchussleJ 
(12OOI1) . 

Latitudes and number of erupting flux tubes 

We assume the emergence rate of flux tubes to be proportional 
to the mean toroidal field given by the dynamo model. 



Bij = B(Ai,tj) 



(19) 
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Latitude (degrees) 



Fig. 3. Number of flux tubes per time interval of 7 days, nj, as a 
function of time. The dashed line shows the value in brackets of 
Eq. (l22l) , and the full line shows nj. 



Fig. 1. Normalised latitudinal profiles of a and the radial shear 




10 20 
A, (degrees) 



Fig. 2. A sample cumulative probability distribution function, 
Qij, as a function of latitude, Ai, at a given time tj. Eor each 
erupting flux tube, a pseudo-random number is chosen in the in- 
terval [0,1]. The latitude Ai at which the flux tube erupts is then 
obtained by the mapping shown by the arrows. 




Fig. 4. Meridional profile of internal rotation used in the stabil- 
ity analysis and flux tube simulations. Contours denote angular 
velocity, normalised to the equatorial value at the surface. 



for a given latitude, Ai, and time, tj, which are discretised into 
time intervals, = tj+i -tj, and latitude intervals, AA = Ai+i -Ai, 
with A = 7i/2 - 6 and -7r/2 < /I < 7r/2. The probability of a flux 
tube at Ai and tj to start rising is assumed to be proportional to 
the local mean field, so that we have 



Bi 



Pij 



(20) 



We denote the number of flux tubes emerging at time tj with 
Hj. To determine the initial latitudes of flux tubes at the bottom 
of the convection zone, nj pseudo-random numbers taken from 
a uniform distribution in the interval [0, 1] are mapped via the 
cumulative probability distribution function (PDE), Qij, based 
on pif. 



Pkj- 



(21) 



The mapping for a sample PDE is illustrated in Eig.|2] The num- 
ber of emerging tubes, uj, is chosen to be proportional to the to- 
tal toroidal flux density at that time, 2/ ^ij, requiring that a total 
number of N tubes emerge per activity cycle, viz. 



^ijBij 



-hO.5 



(22) 



Here the brackets denote the nearest integer to the value of the 
expression they enclose. Eor the case of the Sun, A'^ is set to 2100, 



which constrains the total flux emerging over an activity cycle to 
be of the order of 10^^ Mx, which is comparable with observa- 
tions. Eigure[3]shows nj as a function of time, with B(A, t) taken 
from the solar dynamo model given in Sect. O 

Field strengths of unstable flux tubes 

To determine the field strength ^ft, we consider the stability 
properties of flux tubes (Eerriz-Mas & Schiissler 1995) in the 
middle of the overshoot region (5000 km above the upper bound- 
ary of th e radiative zone), usi ng the non-local mixing length 
model of ISkalev & Stixl (Il99ll) . We consider a rot ation profile 
similar to that obtained using helioseismic data (e.g. JSchou et al.l 
fT998h 



no 



0.9635 ■ 



1 -\- erf 



do 



• (0.0876 cos^ 6 -h 0.0535 cos^ 6 - 0.0182) , (23) 

where Do is the equatorial rotation rate at the surface, and do = 
0.01 SRq is taken as the half-thickness of the tachocline. Lines of 
constant angular velocity are shown in Eig.|4] This 2D profile is 
consistent with the latitude dependence of the radial shear given 
in Eq. ([13]), which we use in the dynamo model. 

Eigure [5^ shows the stability diagram in a plane defined by 
flux tube latitude and field strength. We define Bj(A) to be the 
field strength of unstable flux tubes as a function of latitude. 
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corresponding to a certain growth time, r, of the instability. In 
the course of the ampliflcation of the toroidal magnetic field, in- 
stability sets in when 5ft exceeds a critical value. This value 
also defines the mean field strength in Eq. (O, at which 
the dynamo saturates, taking into account the scaling given by 
Eq. ([TS]) . In most parts of the stability diagram, r decreases with 
increasing field strength. In the aQ. dynamo model, the toroidal 
field amplification occurs on the timescale of the rotational shear 
in the tachocline, which is about a few years. As long as the 
growth time of the flux tube instability is longer than the shear- 
ing timescale, the toroidal field continu es to be amplified withou t 
the flux tube leaving the shear region (llsik & Holzwarthl l2009^. 
As the field strength reaches a level of the order 10^ G, the two 
timescales become comparable and the flux tube is lost from the 
layer. Therefore, we do not consider unstable tubes with growth 
times much longer than the amplification timescale and choose 
the field strengths of the emerging flux tubes given by the curve 
^20o(^), corresponding to r = 200 days. Choosing a growth time 
exceeding this value leads to almost identical results since the 
curves of constant growth time become very close to each other 
near the stability limit (Fig. [5^). We consider the contour line 
^200 (^) of the main region of instability in Fig. [5^, which is a 
single- valued function of latitude. Taking field strengths corre- 
sponding to the low-latitude "island" to the left of the diagram 
does not lead to a significant diff'erence in the emergence lat- 
itudes and tilt angles (Caligari et al. 1995). To avoid the occa- 
sional emergence of flux loops at unrealistically high latitudes 
in the case of the Sun, we assume that unstable flux tubes only 
form if the local mean toroidal field exceeds 400 G (4% of the 
equipartition field strength). 

3.3. The rise of flux loops 

Having determined the times of emergence, initial latitudes, and 
the corresponding field strengths of unstable flux tubes, we sim- 
ulate their rise through the convec tion zone using the code de- 
veloped by lMoreno^ InsertisI (119861) and extended to three dimen- 
sions and spherical geometry bvlCaH sari et al. l (IT99l . In the nu- 
merical scheme, the thin flux tube is described by a string of 
Lagrangian mass elements. 

The initial value for the cross-sectional radius of each flux 
tube is taken to be Rt = 1000 km. For B = 10^ G, this cor- 
responds to a magnetic flux of about 3 x 10^^ Mx, which is 
typical of a bipolar magnetic region (hereafter BMR) of mod- 
erate size at the solar surface. The tube radius is large enough 
for th e rising flux loops to be lar gely unafl^ected by the drag 
force (iD'Silva & ChoudhurillT993l: ICaHgari et al.lll995h . so that 
difl'erential rotation and meridional circulation do not influence 
their motion. Therefore, the dynamics of flux tubes with Rt > 
1000 km is independent of their radius, as long as the thin flux 
tube approximation remains valid. 

The fastest growing modes of unstable flux tubes have az- 
imuthal wavenumbers m = 1 or m = 2. On the other hand, 
the angular separation of large BMRs on the Sun corresponds 
to much larger azimuthal wavenumbers (m = 10 - 60). A possi- 
ble mechanism for limiting the longitudinal separation of BMRs 
is the dynamical disconnection from their magnetic roots at 
depths of a few Mm below the surface soon after emergence 
(ISchiissler & Rempelll2005l) . 

To determine the emergence latitudes and tilt angles, we use 
a table derived from a number of simulations of rising unsta- 
ble flux tubes for various initial latitudes and magnetic field 
strengths. The initial latitudes are in steps of 5° and the corre- 
sponding field strengths are obtained from the curve B2oo(^)- 



Linear interpolation is used to obtain the emergence latitudes 
and tilt angles at intermediate values of initial latitude and 
field strength. Figure |5j) shows a comparison of the dynamo- 
generated toroidal field in the overshoot layer (contour lines) and 
the surface emergence pattern of rising flux loops (dots). The 
overall emergence pattern coincides with the dynamo waves, 
since the poleward deflection of rising flux tubes with field 
strength of the order of 10^ G by the Coriolis force is small for 
solar rotation rates. This result is consistent with the implicit as- 
sumption often made when interpreting solar dynamo models, 
namely that the surface activity pattern closely reflects the dy- 
namo wave pattern. 

4. Surface flux transport 

Times, latitudes, and tilt angles of emerging flux loops obtained 
with the procedures described so far determine the flux input into 
the surface flux transport model, which we consider in this sec- 
tion. The model describes the subsequent evolution of magnetic 
flux at the stellar surface. 

4.1. The surface evolution of magnetic flux 

We follow the evolution of the surface field governed by the 
emergence of BMRs and flux transport by difl'erential rota- 
tion, meridional flow, and (turbulent) difl'usion due to convec- 
tive flows. In contrast to the big flux tubes rising through the 
convection zone, the large-scale surface field described by the 
flux transport model is highly fragmented into small-scale flux 
concentrations. These magnetic elements are vertically oriented 
(owing to buoyancy) and sufficiently small to be passively trans- 
ported by the surface flows via the drag force. 

We therefore assume that the surface field has only a radial 
component, Br, whose evolution as a function of latitude. A, lon- 
gitude, 0, and time, t, i s described by the ni agnetic induction 
equation in the form (cf . iBaumann et al.ll20()4l) 

dBr dBr Id/ \ 

-f = -Q(^)— : + - L(A)BrCOsA] 

at d(p Ri, cos AdA\ J 

^Dhim) + ^ririr) + S(A, 0, 0, (24) 

where 7?^ is the stellar radius, Dh is the term for horizontal dif- 
fusion with uniform turbulent diff'usivity 7]h = 600 km s"\ Dr 
denotes the radial diff'usion term with un iform radial diff'usivity 
rjr = 100 km s"^ (iBaumann et al. 1 1200^ . In addition, S is the 
source term describing the newly emerging BMRs and is the 
surface angular rotation rate as a function of latitude A 

^(A) = 13.38 - 2.30 sin^ ^ - 1.62 sin"^ A degrees day"^ (25) 

( Snodsrass" "1983|). The meridional flo w velocity v 

( Snodgrass & DailevI 119961: iHathawavl Il996h is assumed to 
be 

v(A) = I sin(;r/l//lc) : if W < ^26) 

\ : otherwise, ^ ' 

wher e f;o = 1 1 m s~^ a nd Ac = ±75° (cf. Ivan Ballegooijen et al.l 
ll998l:lBaumann et al.ll2004l) . 

The numerical solution of Eq. (l24l) is carried out by repre- 
senting the magnetic field as an expansion in ter ms of spherical 
harmonics with a maximum degree of ^ = 63 (iBaumann et al.l 
|2004|) . This corresponds approximately to the observed size of 
the supergranules (-30 Mm) on the Sun. 
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Fig. 5. (a) Stability diagram for thin magnetic flux tubes in the middle of the solar overshoot region, for Sun-like diflferential rotation 
(given by Eq. [23l). White regions represent stable flux tubes, while shaded areas indicate instability. Eor light- shaded areas, the 
fastest growing mode has azimuthal wave number m = 1, while for dark-shaded areas it is m = 2. Contours show lines of constant 
growth times (e-folding times), labelled in units of days, (b) Time-latitude diagram of the dynamo-generated mean toroidal magnetic 
field B(A, t) at the bottom of the convection zone (contours at ±0.04, ±0.1, ±0.3, ±0.5, ±0.7, ±0.9 of the maximum value of B, solid 
for positive and dashed for negative values) and of the flux loops emerging at the surface (dots) for a Sun-like star, (c) Time-latitude 
diagram of the longitudinally averaged radial surface magnetic field for the solar model, (d) Time variation of the total unsigned 
surface magnetic flux. The values are averaged over 27 -day time intervals. 
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Fig. 6. Time variation of the emerging magnetic flux integrated 
over the surface, and summed over 27-day time intervals. The 
corresponding average unsigned field strength is indicated on the 
right-hand vertical axis. 



4.2. Treatment of source BMRs 

The flux density of a newly emerged BMR is written in the form 

Br(A,cp) = B\A,cp)-B-(A,(P). (27) 

The emergence longitu des of the BMRs ar e ran domly 
distr ibuted. Fo l lowin g Ivan Ballegooijen et al.l (Il998l) and 
.Baumann et a D ll2004 . we assume that the unsigned field 
strength of the two polarities is given by 



B-(A, (p) = BoQxp 



St 



(28) 



where I3+(A, 0) are the heliocentric angles between a given posi- 
tion (A, (p) and the centres of the positive and negative polarities, 
(A+, (p+). The initial characteristic angular width of each polarity, 
^in, and the angular separation of the centres of the two poles, Aj3, 
are related by = OAAp. The max imum field strength is set 
to 250 G (see lBaumann et al.|[2QQ4l) . 

The BMRs, which are smaller than the grid cell size (~ 1 
square degree) of the simulations, are taken into account by con- 
sidering these regions only at a later stage in their evolution, after 
having expanded (by turbulent diff'usion) to an angular width of 
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^0 = 4°. The description of the two polarities for such a BMR is 
thus given by 



Bo\ 1 1 exp 



^1 



(29) 



In our numerical treatment, the spherical harmonic coefficients 
defining the field distribution of the newly emerged BMRs are 
multiplied by a spatial filter of the form exp[ - + l)/4j 
(Ivan Ballegooijen et al.ll 19981) . in order to diminish the effect of 
ringing (Gibbs phenomenon) caused by the truncation of the ex- 
pansion in spherical harmonics. 

4.3. Size distribution and the emerging magnetic flux 

The total size of a newly emerged BMR is related to the helio- 
centric angular separation of the centres of the two poles, AjS, 
which, in our model, ranges from 3.5° to 10° in steps of 0.1°. 
The corresponding range of areas is 30-250 square degrees, and 
the resulting range of magnetic ffux is 2.3 x 10^^ - 3.5 x 10^^ 
Mx. The areas A of the bipolar source regions are determined 
randomly with a numbe r distribution derived f rom solar obser- 
vations, N(A) ~ ( Sc hriiver & Harvevl[l994l) . 

For the solar model, the time series of emerging magnetic 
ffux, summed over 27 days (one full solar rotation), is given in 
Fig.[6j It shows a cyclic variation with a period of about 1 1 years, 
closely following the underlying oscillatory dynamo model. 

The evolution of the surface ffux is given in a supplementary 
animation (see the online appendix, Fig. lA.U top panels), which 
shows the weakening and cancellation of the polar ffeld, and the 
gradual build-up of the polar ffeld of the following cycle. The 
time-latitude diagram of the longitudinal average of the signed 
ffux density is shown in Fig. [5]:. The extreme v alues of the un- 
signed ffux density are lower than those found by lBaumann et al.l 
(|2004|) for a similar model. The reason is that the tilt angles re- 
sulting from the simulated rise of ffux tubes ar e systematically 
smalle r than the latitude dependence assumed by lBaumann et al.l 
(|2004|). but actually pro vide a better match to the observations 
(see lCaligari et al.iri995l Fig. 12). The corresponding time vari- 
ation in the total unsigned ffux (and the equivalent average un- 
signed ffeld strength) is shown in Fig. [5]l. The activity minima 
are less marked than those of the emerging ffux (Fig. |6]l; this 
is due to the phase difference of about half a period between the 
variations in the polar ffux and the ffux emerging at low latitudes. 



5. Rapidly rotating Sun-lil<e stars 

As a ffrst non- solar application of the combined model of mag- 
netic ffux generation and transport, we consider rapidly rotating 
stars of solar internal structure. 

For simplicity (and due to lack of information), we assume 
the same values for the turbulent magnetic diffusivities for the 
dynamo region, the convection zone, and the surface as in the 
solar model discussed in the previous sections. We also use the 
same proffles for the rotational shear, i.e., we assume that the dif- 
ference between the minimum and maximum angular velocities 
in the convective envelope. All, is independent of the equatorial 
rotation rate at the surface, Q^. This is motivated by observa- 
tions of cool stars indicating a rather weak dependen ce of AQ at 
the st ellar surface on (e.g., Barnes et al. 2005; Do nahue et al.l 
Il996l) . This assumption corresponds to a decrease in AQ/Q* 
by a factor ^lo/Il* with respect to the solar value, so that the 



constants on the right hand side of Eq. (1231) have to be multi- 
plied by Qq/II^. Because AD. is kept constant, the surface ffux 
transport equation (Eq.|24l), the rotation proffle (Eq. |25}, and Rq 
(Eq. [9]) remain unchanged. We also adopt a solar- type internal 
rotation proffle in the present section, owing to the lack of obser- 
vational information on stellar internal rotation proffles. When 
results from asteroseismology become available, they can be di- 
rectly incorporated into the model. 

In the dynamo mod el, we assume that ao scales with 
(iKrause & RadleJl980l) , so that Rq. oc in Eq. ((8]), because rj is 
not changed. The dynamo cycle period, Pcyc, then decreases with 
increasing Q^. We note, however, that this scaling might not be 
valid for the very most rapid rotators since the a-effect driven by 
helical turbulence saturate s when the cyclonic rotation of rising 
fluid parcels exceeds n/2 (iParkejl 19821: iRiidiger & Kitchatinovl 
119931) . 

The total number of erupting flux tubes per activity cycle, N, 
is also scaled with Q^, in accordance with the observed relation 
between the rotation rat e and the level of magnetic activity (e.g., 
iMontesinos et al.ll200lh . 



5.1. The case ofProt = 9 d 

We ffrst consider a Sun-like star with an equatorial rotation rate 
of Di, = 2.1 Dq, corresponding to a period of about 9 d. Figure|7t 
shows the stability diagram of ffux tubes located in the middle 
of the overshoot region. 

The critical ffeld strength for the Parker instability is larger 
than in the solar case, since rotation has a stabilising effect on 
the ffux tubes. With increasing D^,, angular momentum conser- 
vation suppresses the tube displacements perpendicular to the 
rotation axis, so that the magnetic buoyancy needs to be stronger 
than in a similar situation in the Sun. 

In determining the function 520o(^), we excluded the 'is- 
lands' of instability seen in Fig.|7^. Test simulations showed that 
the emergence latitudes and tilt angles of rising ffux tubes do not 
differ signiff cantly between the islands and the neighbouring part 
of the main region of instability. 

The time-latitude diagram of the dynamo-generated toroidal 
ffeld at the bottom of the convection zone and the locations of 
ffux emergence are combined in Fig.|7J). They show a cyclic vari- 
ation with the dynamo period of about ffve years. The poleward 
deff ection of rising ffux tubes is slightly stronger than in the case 
of solar rotation. The effect responsible is the component of the 
Coriolis acceleration directed towards the rotation axis (resulting 
from a retrograde ffow along the ffux tube required by angular 
momentum conservation), which is proportional to the rotation 
rate. 

The emerging magnetic ffux, summed over one-month inter- 
vals is shown in Fig. (8] While the total ffux emerging per cycle 
is about ffve times higher than in the solar case, the relative vari- 
ation between the minima and maxima of the emerging ffux is 
considerably weaker. This is due to the larger overlap between 
consecutive cycles, which is a consequence of stronger dynamo 
excitation. Figure |7]: shows that also the polar ffeld maxima are 
higher than in the case of solar rotation (note the different sat- 
uration level of the colour table from Fig. |5]:). This is due to 
1) a higher ffux emergence rate and 2) the bigger tilt angles of 
the emerging BMRs because of the stronger Coriolis force expe- 
rienced by the horizontally expanding crest of the rising tubes. 
The animation of the surface ffeld distribution is given in a sup- 
plementary animation (see the online appendix, Fig. lA.U middle 
panels). Figure |7]l shows that the cycle signal introduced by the 
underlying dynamo is largely washed out in the time evolution 
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Fig. 7. Same as Fig.O but for a Sun-like star with Prot = 9 d. 



Magnetic flux emerging per month 




latitude than their starting location (ISchtissler & Solankilll992h . 
1 2Q The initial latitudes are in the range of roughly 0°-25°, while the 
^ emergence latitudes are about 23°-40°, so that there is no flux 
emergence in an 'avoided zone' up to ±23° latitude from the 
§ equator. Figure |9j) shows the dynamo-generated toroidal mag- 
^ netic field at the bottom of the convection zone together with the 
distribution of flux tubes emerging at the surface. It is evident 
that the poleward deflection of flux tubes causes a significant de- 
viation between the dynamo waves in the interior and the pattern 
of surface emergence. 



Fig. 8. Same as Fig.[6l but for a Sun-like star with P^ot = 9 d. 



of the total magnetic flux at the surface. This is caused by 1) a 
larger overlap between consecutive cycles owing to a stronger 
a-eflfect, and 2) strong polar magnetic fields in antiphase with 
the cycle of emerging flux. 

5.2. The case ofProt = 2 d 

For a star with a rotation rate of = ISIl©, corresponding 
to a period of two days, the stabilising eff'ect of rotation on the 
flux tubes is even stronger than in the previous case. The stabil- 
ity diagram given in Fig.[9t shows that the critical field strength 
increased to a few times 10^ G. As before, we exclude the in- 
stability islands in determining the field strengths of rising flux 
tubes. 

For such a large rotation rate, the flux tubes rise almost par- 
allel to the rotation axis and emerge at a significantly higher 



For such a rapidly rotating star, the tilt angles of emerging 
bipolar regions are much larger (around 35°) than in the solar 
case. The relatively narrow latitudinal extension of flux emer- 
gence and the large tilt angles lead to a pair of almost non- 
migrating latitudinal belts of opposite magnetic polarity in each 
hemisphere (see Fig.[9j: and the supplementary animation in the 
online appendix. Fig. lA.U bottom panels). The net poleward 
transport of flux from the higher-latitude bands produces polar 
fields reaching field strengths of up to 35 G. Compared to the 
Prot = 9 d case, there is smaller overlap between the emergence 
patterns of consecutive cycles. Consequently, the correspond- 
ing surface-integrated unsigned magnetic flux (Fig. |9ji) shows 
a clearer cycle (mainly determined by the activity bands), with a 
period of about two years, an amplitude of about 1 .5 x 10^^ Mx, 
and a minimum-to-maximum flux variation by a factor of about 
two. 
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Fig. 9. Same as Fig.O but for a Sun-like star with Prot = 2 d. In panel (b), the dots are shown at one-month intervals for clarity. 



6. Rapidly rotating K stars 

6.1. KO dwarf 

We consider a KOV star at the zero-age main sequence, with a 
rotation period of 2 days. Such a star has a smaller radius and a 
relatively deeper convection zone than the Sun. We use a stellar 
model with R^, = O.lSS Rc?) and the bottom of the co nvection 
zone located at 0.6777?^ (iHolzwarth & Schusslejl200lh . 

The Reynolds numbers for differential rotation, Rq, and the 
Qf-effect, Ra, have to be adapted to the new value of the length 
scale, namely the radial distance of the convection zone lower 
boundary from the center. Scaling with D^/Dq, we obtain ao = 
152 cm s"^ The resulting dynamo number R^Rq is about a factor 
of four s malle r than in the case of the Sun-like star with P^ot = 
2 d (Sect. l5.2l) , which leads to a somewhat longer magnetic cycle 
period of about 6.5 years. The number of emerging flux tubes per 
cycle is set to the same value as in the Sun-like star with the same 
rotation period. The emergence latitudes and tilt angles of the 
BMRs are determined with the same procedure as in the previous 
sections, considering the stability properties of flux tubes in the 
mid-overshoot region of the KOV star (Fig.fTOh). 

The comparison of the dynamo waves of the deep-seated 
toroidal field with the emergence pattern at the surface is shown 
in Fig.fTOb. The diflTerence between the emergence latitudes and 
the initial latitudes of flux tubes is larger by about 10° than for 
the solar- type star with the same rotation rate. This is due to a 
simple geometric eflTect: the KOV star has a deeper convection 
zone (in terms of fractional radius), so that a flux loop rising 



parallel to the rotation axis emerges at a higher latitude than in 
the case of a shallower convection zone. 



The time-latitude diagram of the longitudinally averaged 
magnetic field and the time variation of the integrated unsigned 
surface flux are shown in Fig. [TOb . Owing to the rapid rotation, 
the emerging BMRs have large tilt angles between 40° and 55°. 
This leads to pronounced branches of follower polarity drift- 
ing poleward via difl'usion and meridional circulation. The cycle 
signal is similar in strength to the corresponding Sun-like case 
(Sect. 15.2b , again due to the small cycle overlap of the emer- 
gence patterns and the large tilt angles leading to the formation 
of activity bands (Fig.fTOb). 



Three snapshots of the flux distribution at the stellar surface 
are shown in Fig. [TT] Near an activity minimum (left panel), 
mid-latitude BMRs of the newly started cycle lead to an un- 
even longitudinal distribution of magnetic flux, owing to the rel- 
atively low emergence frequency. In the maximum phase (mid- 
dle panel), highly tilted BMRs lead to the formation of a pair 
of latitudinal belts of opposite polarity in each hemisphere. The 
higher-latitude belt provides the source for the respective po- 
lar field of the same polarity. After seven months (right panel), 
flux emergence has decreased in parallel to the evolution of the 
deep-seated toroidal field generated by the dynamo, while the 
polar field has reached its maximum strength. As a result, the 
surface-integrated unsigned magnetic flux exhibits a cyclic vari- 
ation with a period of about 2.75 yr (Fig.fTOh). 
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Fig. 10. Same as Fig.S but for a KOV star with P,ot = 2d. 



6.2. K1 subgiant 

As a further example, we consider a Kl subgiant rotating with 
a period of 2.8 days. A star of this type is the magnetically ac- 
tive component of HR 1099, the most thoroughly studied RS 
CVn-type binary system. In t he structure model adopted here 
(iHolzwarth & Schussle^l200lb . the stellar mass is l.SM©, the 
age is about 2.4 Gyr, the stellar radius is 7?^ = 3.447?q, and the 
bottom of the convection zone is at 0.357?^. These parameters 
still allow flux tubes to rise up to the surface. If the bottom of the 
convection zone were located deeper than about 0.307?*, buoy- 
ant flux tubes entering the convection zone would be trapped 
within the star as a consequen ce of the magnetic curvature force 
(IHolzwarth & Schussledl200l1) . 

We assumed rigid rotation in the calculations regarding the 
instability and the rise of magnetic flux tubes, and the sur- 
face flux transport. This assumption is sup ported by the very 
weak su rface shear observed in HR 1099 (iDonati et al.l 120031: 
iPetit et al. 2004). Furthermore, the eflfects on flux tube dynam- 
ics of difl'erentially rotating convection zones in post-main se- 
quence stars is yet to be explor ed, as an extension of the work by 
IHolzwarth & Schusslej (l200ll) . The linear stability diagram for 
flux tubes in the middle of the overshoot region (about 0.337?*) 
is shown in Fig.[T2h. 

For the dynamo model, Rq has been rescaled using the radius 
ro = IARq of the star's radiative core and the surface latitudinal 
shear of about 0.87°d"\ fol lowing the measured value with the 
lowest variance, as given by lPetit et al.l (l2004l) . The resulting dy- 
namo number leads to an activity cycle with a period of about 4 




Fig. 11. Snapshots of the magnetic field distribution at the sur- 
face of the KOV star with Prot = 2 d. The colour table shows the 
field strength with a saturation level at ±500 G. Left: the activity 
minimum phase. Middle: the activity maximum phase. Right: 
the beginning of the decay phase following the activity maxi- 
mum, when the polar field reaches its maximum strength. The 
inclination of the rotation axis with respect to the line of sight is 
30°. 



yr, compared to the ob served cycle periods of 5.3 yr and 16 yr 
(iBerdvugina & Henrvll2007l) . 

Figure [T2b shows that the dynamo waves and the emergence 
pattern have no spatial correlation. As a result of the d eep c on- 
vection zone and the geometric eflfect described in Sect. 16.11 the 
average latitude of emergence is shifted to even higher latitudes 
than in the case of main- sequence stars. The emergence belt cen- 
tred around A ^ 62° has a latitudinal width of only about 7 de- 
grees. This is caused by the projection of a low-latitude angu- 
lar range at the bottom of the convection zone, onto the surface 
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and along the direction parallel to the rotation axis, producing a 
much narrower range at high-latitudes. 

Emerging BMRs have tilt angles between 5° and 8°, which 
are significantly smaller than those for the rapidly rotating main- 
sequence stars considered in the previous sections. This can be 
explained in the present case by the rising flux tubes spending a 
much longer time (about five years) in the convection zone, so 
that the magnetic tension force can more efl&ciently act against 
the twisting eflTect of the Coriolis force on the horizontally ex- 
panding loop. 

The time-latitude diagram of the longitudinally averaged 
magnetic field strength (Eig.[T2b) is similar to that of the rapidly 
rotating Sun-like star, apart from a shift of the patterns to higher 
latitudes. The total unsigned surface flux variation (Eig. \V2h ) 
shows a rather weak cyclic signal upon which irregular varia- 
tions are superimposed. The narrow latitude range of flux emer- 
gence, as well as the weak transport eflfects by surface flows and 
diflTusion in this region, amplify the eff'ect of the random compo- 
nent of the surface flux injection process. This results in rather 
large fluctuations of the total unsigned surface flux. 

Eigure [T2b also indicates a rather weak field in the polar 
caps above -75° latitude, although polarity reversals still oc- 
cur. This weakness is due to the relatively small tilt angles of the 
BMRs and the long diflfusion timescale (proportional to R^)\ the 
meridional flow vanishes beyond ±75°, so that further poleward 
flux transport occurs only by means of the relatively ineflftcient 
diflfusion process. Inefl&cient diflfusion and small tilt angles lead 
to an intermingling of polarities in the narrow flux emergence 
belts, as shown in Eig. [131 Therefore, the resulting poloidal field 




Fig. 13. Eield strength distribution at an activity minimum (left 
panel), and at a maximum (right panel), for the Kl subgiant star. 
The inclination of the rotation axis with respect to the line of 
sight is 30°. The colour table shows the field strength with a 
saturation level at ±300 G. 

produced by the diffusing BMRs is weak and far from a simple 
dipole configuration. 

7. Discussion 

The combined model of magnetic field generation and transport 
provi des two major improvements compared to pr evious models 
(e.g. jBaumann et al ] |2004I: iHolzwarth et aDl2006h : 

1. The cycle properties are determined by an underlying dy- 
namo model, whereas in previous studies either the butterfly 
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diagrams were prescribed arbitrarily, or observed emergence 
patterns were used. We thus establish the nontrivial link be- 
tween the deep-seated field generation and surface flux trans- 
port. 

2. The tilt angles and emergence latitudes of bipolar magnetic 
regions at the surface are not prescribed, but consistently 
determined from simulations of rising flux tubes. This may 
strongly aflTect the surface distribution of magnetic flux, and 
particularly the strength of the polar magnetic fields. 

In this first exploratory study, our main motivation was to 
evaluate the eff'ects of combining models for the dynamo, flux 
tube stability and dynamics, and surface transport, rather than to 
constrain dynamo models for cool stars. We therefore adopted a 
rather simple thin-layer dynamo model with difl'erential rotation 
in radius, which does not include the probably important efl'ects 
of latitudinal difi^erential ro tation and meridional circulation. We 
note also that lSpruil (1201 Ol) argues that the radial angular veloc- 
ity gradient in the solar tachocline is not capable of building up a 
sizeable toroidal field. Our model can certainly be utilised to test 
more elaborate dynamo models, such as Bab cock-Leighton-type 
flux transport dynamos (ICharbonneaull2Q05l) . 

For solar- type and KO main- sequence stars, we find that 
the polar magnetic fields for rotation periods of 10 d and 2 d 
are stronger than in the solar case, indicating the impo rtance 
of surf ace transport eff'ects ( Schrijver & Title 2001; Is ik et al.l 
I2007bh : the emergence and meridional transport of BMRs at 
mid-latitudes with tilt angles as high as 35° lead to strong (in 
comparison to the Sun) unipolar flux regions around the poles. 

For the solar- type model with Pj-ot = 10 d, the surface flux 
transport blurs the periodic signal of the oscillatory dynamo in 
the surface-integrated unsigned flux. This suggests that for rapid 
rotators, a cyclic dynamo may actually underly an apparently ir- 
regular variation in activity indicators. This could possibly be 
the case f or HP 190 7 71, a Sun-like star with P^ot = 8.8 d ob- 
served by iPetit et all (l2009l) . Using the Doppler-imaging tech- 
nique, they report that the total magnetic energy at the surface 
does not change over the observed period, although global polar- 
ity reversals in the radial and azimuthal field components occur. 
Although considering only a radial field at the surface, our mod- 
els show that a roughly constant total magnetic flux at the surface 
may arise from the combined eff'ects of dynamo cycle overlap, 
large tilt angles, high flux emergence frequency, and meridional 
transport. 

For all models with P^ot = 2 d, we found that the dynamo 
wave pattern and the emerging surface flux are completely dif- 
ferent from each other, owing to the strong poleward deflection 
of rising flux loops by the Coriolis force. Consequently, for rapid 
rotators the observed activity patterns on the surface may not 
necessarily be taken to represent the spatio-temporal distribution 
of magnetic field in the deep interior. 

The results presented here are certainly rather preliminary, 
owing to the very simple dynamo model used and the insuf- 
ficient information we have at present about various important 
parameters: 

- The dependence on rotation rate and stellar structure of a- 
eflfect, internal rotation profile, and the meridional flow ve- 
locity are largely unknown. 

- Various complications regarding the formation of flux tubes, 
the growth of the buoyancy instability, and the rise through 
the convection zone are omitted, because the connection be- 
tween dynamo-gener ated field and flux tube dynami cs is not 
well understood (e.g. JSchiissler & Ferriz-Mas1l2003l) . 



- Reliable models for the structure of convective overshoot re- 
gions do not exist. 

- It is unclear whether the solar value of rjh = 600 km^ s"^ 
for the turbulent surface diff'usivity is valid for stars of dif- 
ferent structure, rotation, and fraction of surface covered by 
magnetic flux. 

For the solar- type stars with P^ot = 10 d and 2 d and the 
KOV star with Prot = 2 d, low/mid-latitude flux emergence and 
strong polar magnetic fields co-exist. However, the lack of mag- 
netic regions for \A\ < 60° for the Kl subgiant is in disagree- 
ment with obser vations for which active regions are detecte d 
at lower latitudes (IVogt et al.lll999l:IStrassmeier & Bartusll2000h . 
This could be due to an underestimation of the required field 
strength for flux tube instability and rise: if the subadiabaticity 
in the overshoot region of the subgiant were higher than given 
by the stellar model used, unstable rising flux tubes would have 
a stronger field and thus be more buoyant, leading to less pole- 
ward deflection by the Coriolis force and their emergence at 
lower latitudes. Another possible explanation is non-solar merid- 
ional flow profiles (e.g., with large multiple cell structures), 
which can have a large impact on the cycle pro perties, par- 
ticula rly for advection-dominated dynamo models (iJouve et al.l 
I2OIOI) . Furthermore, these patterns could also provide equator- 
ward flux transport. 

Our assumption that the number of flux tubes emerging 
per activity cycle is proportional to the rotation rate is qualita- 
tively consistent with our scaling of the a-efl'ect with rotation 
rate. However, there are energetic and geometrical constraints 
on the maximum amount of magnetic flux that can be generated 
in rapid-rotato r convection zones with the known mechanisms 
(lRempelll200'8b . 



8. Conclusions 

We have presented a combined model including 1) a thin-layer 
dynamo with an a-eff'ect due to magnetic buoyancy instability 
and radial shear according to helioseismic results, 2) the stability 
analysis and dynamical simulation of flux tubes rising through 
the convection zone, and 3) surface flux transport by diff'erential 
rotation, meridional flow, and turbulent diflfusion that reproduce 
the basic observed features of the solar cycle (such as the butter- 
fly diagram of flux emergence, the evolution of the large-scale 
surface flux distribution, and the polar field reversals). 

Extending the model to rapidly rotating Sun-like and K-type 
main- sequence stars, we have found that the latitude distribution 
of emerging flux at the surface may diff'er significantly from that 
of the deep-seated toroidal field generated by the dynamo pro- 
cess. A strong overlap of the cycles may hide the cyclic signal 
in activity indicators related to the integrated total surface flux. 
Low- to mid-latitude flux emergence and activity can co-exist 
with strong polar fields ('polar spots') arising from poleward sur- 
face flux transport. 

In the case of a Kl subgiant, we have found that the large 
depth of the convection zone and long diff'usion time lead to 
belts of intermingled mixed-polarity field at high latitudes and 
relatively weak polar field. 

For improved realism of the model, a two-dimensional flux- 
transport dynamo model will be incorporated at our next pro- 
posed developed stage. Information about the stellar interior and 
the surface diff'erential rotation, meridional flow, and internal 
structure can be included in the model as it becomes available 
by means of asteroseismology. 
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Fig. A.l. The surface distributions of magnetic flux for the Sun- 
Hke star with Prot = 27 d (top panels), Prot = 9 d (middle panels), 
and Prot = 2 d (bottom panels), near the activity minimum (left 
panels) and the maximum phases (right panels). The colour scale 
for the magnetic field strength saturates at ±5, 30, and 75 G for 
top, middle, and bottom panels, respectively. The correspond- 
ing time-latitude diagrams for the azimuthally averaged mag- 
netic field strength are shown in Figs. [5]:, [7]:, and|9]:. 



Appendix A: Animations 

We provide three animated GIF files showing the simulated evo- 
lution of surface magnetic flux for three rotation rates of the Sun- 
like star, available on-line. 

The animation Bsurf _27d . gif shows the evolution of mag- 
netic field at the surface for the Sun-like mod el (Se ct. 14.31 ). Still 
frames from the animation are shown in Fig. lA.lK top panels). 
The simulation begins at a cycle phase close to the polar field 
maximum with negative polarity at the north pole. The overlap- 
ping of two activity cycles is visible, i.e., high-latitude BMRs of 
the next cycle emerge, while BMRs of the ending cycle continue 
to emerge at low latitudes. In the course of the cycle, the polar 
field is weakened by cancellation with opposite-polarity (nega- 
tive) flux transported from the active lower latitudes and even- 



tually reverses around activity maximum. Afterwards, BMRs of 
the new cycle gradually build up positive flux at the north pole. 

The animation Bsurf_9d.gif shows the variation in mag- 
netic field at the surface of the Sun-like star with Prot = 9 d 
(Sect. l5.ll) . Still frames from the animation are shown in Fig. lA.ll 
(middle panels). Although the flux emerges within a rather re- 
stricted latitude range, it is spread over the entire stellar surface 
by diff'usion and the large-scale flows. There are unipolar mag- 
netic regions near the poles and mixed-polarity regions at low to 
middle latitudes. 

The animation Bsurf_2d. gif shows the evolution of mag- 
netic field at the surface of the Sun-like star with Pj-ot = 2 d 
(Sect.[5]2]). Still frames from the animation are shown in Fig. lA.ll 
(bottom panels). The latitude range of flux emergence is smaller 
than in the previous cases, and the tilt angle is significantly 
larger. As a result, two bands of magnetic regions of opposite 
polarity form, predominantly during the activity maximum. 



